JLAB-THY-03-03, ADP-03-107/T545 



Bound Nucleon Form Factors, Quark-Hadron Duality, and 

Nuclear EMC Effect 

K. Tsushima 1 *, D.H. Lu 2 , W. Melnitchouk 3 , K. Saito 4 , and A.W. Thomas 5 

1 Department of Physics and Astronomy, University of Georgia, Athens, Georgia 30602, USA 
2 Department of Physics and Zhejiang Institute of Modern Physics, Zhejiang University, 

Hangzhou 310027, China 
3 Jefferson Lab, 12000 Jefferson Avenue, Newport News, VA 23606, USA 
^Physics Division, Tohoku College of Pharmacy, Sendai 981-8558, Japan 
^Department of Physics and Mathematical Physics, and Special Research Centre for the 
• Subatomic Structure of Matter, Adelaide University, SA 5005, Australia 

(N 
£ 

We discuss the electromagnetic form factors, axial form factors, and struc- 
ture functions of a bound nucleon in the quark-meson coupling (QMC) model. 
Free space nucleon form factors are calculated using the improved cloudy bag 
model (ICBM) . After describing finite nuclei and nuclear matter in the quark- 
based QMC model, we compute the in-medium modification of the bound 
£ — . \ nucleon form factors in the same framework. Finally, limits on the medium 

modification of the bound nucleon F<i structure function are obtained using 
«) ■ the calculated in-medium electromagnetic form factors and local quark-hadron 

■ duality. 
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I. INTRODUCTION 



O . 

Partial restoration of chiral symmetry in a nuclear medium, to which the reduction in 
the mass of a bound nucleon is sometimes ascribed, plays a key role in understanding the 
medium modification of bound nucleon (hadron) properties. At relatively high energies 
and/or temperature and/or densities, quark and gluon degrees of freedom are expected to 
be efficient in describing physical phenomena according to perturbative QCD. However, it 
is not at all obvious whether such degrees of freedom are indeed necessary or efficient in 
describing low energy nuclear phenomena, such as the static properties of finite nuclei. In 
this article, we demonstrate that the quark degrees of freedom do indeed seem to be necessary 
to understand recent polarization transfer measurements in the 4 He(e, e'p) 3 H reaction [1,2], 
which cannot be explained within the best existing treatments of traditional physics (solely 
based on hadronic degrees of freedom). 

Over the past few years there has been considerable interest in possible changes in bound 
nucleon properties in a nuclear medium. There is a significant constraint on the possible 
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change in the radius of a bound nucleon based on y-scaling — especially in 3 He [3] . On the 
other hand, the space (time) component of the effective one-body axial coupling constant is 
known to be quenched [4] (enhanced [5]) in Gamow- Teller (first-forbidden) nuclear (3 decay, 
and a change in the charge radius of a bound proton provides a natural suppression of the 
Coulomb sum rule [6] . One of the most famous nuclear medium effects — the nuclear EMC 
effect [7] , or the change in the inclusive deep-inelastic structure function of a nucleus relative 
to that of a free nucleon — has stimulated theoretical and experimental efforts for almost 
two decades now which seek to understand the dynamics responsible for the change in the 
quark-gluon structure of the nucleon in the nuclear environment [8]. 

Recently the search for evidence for some modification of nucleon properties in medium 
has been extended to electromagnetic form factors, in polarized (e, e'p) scattering exper- 
iments on 16 O [9] and 4 He [1,2]. These experiments measured the ratio of transverse to 
longitudinal polarization of the ejected protons, which for a free nucleon is proportional to 
the ratio of electric to magnetic elastic form factors [10], 



G\ _ P x E e + E' e 
G P m " P' z 2M N 



tan(0 e /2) . (1) 



Here P' x and P' z are the transverse and longitudinal polarization transfer observables, E e 
and E' e the incident and recoil electron energies, 9 e the electron scattering angle, and M N 
the nucleon mass. Compared with the traditional cross section measurements, polarization 
transfer experiments provide more sensitive tests of dynamics, especially of any in-medium 
changes in the form factor ratios. The feasibility of this technique was first demonstrated 
in the commissioning experiment at Jefferson Lab on 16 O [9] at Q 2 = 0.8 GeV 2 . In the 
subsequent experiment at MAMI on 4 He [1] at Q 2 0.4 GeV 2 , and at Jefferson Lab at 
Q 2 = 0.5, 1.0, 1.6 and 2.6 GeV 2 , which had much higher statistics, the polarization ratio in 
4 He was found to differ by ~ 10% from that in 1 H. 

Conventional models using free nucleon form factors and the best phenomenologically 
determined optical potentials and bound state wave functions, as well as relativistic cor- 
rections, meson exchange currents, isobar contributions and final state interactions [11-14], 
fail to account for the observed effect in 4 He [1,2]. Indeed, full agreement with the data 
was only obtained when, in addition to these standard nuclear corrections, a small change 
in the structure of the bound nucleon, which had been estimated within the quark-meson 
coupling (QMC) model [15-21], was taken into account. The final analysis [2] seems to favor 
this scenario even more, although the error bars may still be too large to draw a definite 
conclusion. 

On the other hand, there has recently been considerable interest in the interplay between 
form factors and structure functions in the context of quark-hadron duality. As observed 
originally by Bloom and Gilman [22] , the F 2 structure function measured in inclusive lepton 
scattering at low W (where W is the mass of the hadronic final state) generally follows a 
global scaling curve which describes high W data, to which the resonance structure function 
averages. Furthermore, the equivalence of the averaged resonance and scaling structure 
functions appears to hold for each resonance region, over restricted intervals of W, so that 
the resonance-scaling duality also exists locally. These findings were dramatically confirmed 
in recent high-precision measurements of the proton and deuteron F 2 structure function at 
Jefferson Lab [23,24], which demonstrated that local duality works remarkably well for each 
of the low-lying resonances, including surprisingly the elastic, to rather low values of Q 2 . 
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In this article we first briefly review how finite nuclei and nuclear matter are treated in the 
quark-based QMC model [20,21]. We then discuss the modification of the electromagnetic 
and axial form factors of a bound nucleon in the same model. Finally, using the concept 
of quark-hadron duality and the calculated bound nucleon electromagnetic form factors, 
we extract the F 2 structure function of the bound nucleon [25]. To the extent that local 
duality is a good approximation, the relations among the nucleon form factors and structure 
functions are model independent, and can in fact be used to test the self-consistency of the 
models. We find that the recent form factor data for a proton bound in 4 He [1,2] place strong 
constraints on the medium modification of inclusive structure functions at large Bjorken-x. 
In particular, they appear to disfavor models in which the bulk of the nuclear EMC effect 
is attributed to deformation of the intrinsic nucleon structure off-shell - see e.g. Ref. [26]. 

This article is organized as follows. In Section II we briefly review the treatment of 
finite nuclei in the QMC model [20,21]. Then, in Section III, we discuss the in- medium 
modification of bound nucleon electromagnetic form factors in the QMC model [15-18] as 
inferred from the recent polarization transfer experiments, as well as that of the axial form 
factor Ga{Q 2 ) [27]. In Section IV quark-hadron duality is used to relate the observed 
form factor modification to that which would be expected in the deep-inelastic F 2 structure 
function [25]. Finally, we summarize our findings in Section V. 



II. FINITE NUCLEI AND NUCLEAR MATTER IN THE QMC MODEL 

In this Section we briefly review the treatment of finite nuclei and symmetric nuclear 
matter in the QMC model [19-21]. We consider static, spherically symmetric nuclei, and 
adopt the Hartree, mean-field approximation, ignoring the pNN tensor coupling as usually 
done in the Hartree treatment of quantum hadrodynamics (QHD) [28] (see Refs. [20,21] for 
discussions about the pNN tensor coupling). 

Using the Born-Oppenheimer approximation, mean-field equations of motion are derived 
for a nucleus in which the quasi-particles moving in single-particle orbits are three-quark 
clusters with the quantum numbers of a nucleon. A relativistic Lagrangian density at the 
hadronic level can then be constructed [20,21], similar to that obtained in QHD [28], which 
produces the same equations of motion when expanded to the same order in velocity: 



r N e 
3 i /i i ~N\ 



i-y-d- M* N {a) ~ ( + 9 P fb(r) + -(1 + r 3 N )A(f) ) 7o 



4[(V<7(rO) 2 + mlam + i[(Vc(f)) 2 + m>(rf] 

+ \[{Vb{^f + mlb{rf} + \{VA{r))\ (2) 

where ^jv(r) and b(f) are respectively the nucleon and p meson (the time component in the 
third direction of isospin) fields, while m a , m u and m p are the masses of the a, u and p 
meson fields. g u and g p are the u-N and p-N coupling constants which are related to the 
corresponding (u, d)-quark-u;, g%, and (u, d)-quark-p, g q , coupling constants by g w = 3g^ 
and g p = g q [20,21]. (Hereafter, we will denote the light quark flavors by q = u, d.) The field 
dependent a-N coupling strength predicted by the QMC model, g a (c), which is related to 
the Lagrangian density of Eq. (2) at the hadronic level, is defined by: 
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M* N (a) = M N - g a (a)a(f) . (3) 

Note that the dependence of these coupling strengths on the applied scalar field must be 
calculated self-consistently at the quark level [20,21]. From the Lagrangian density in Eq. (2), 
a set of equations of motion for the nuclear system can be obtained: 

[< 7 • d - M* N {<j) - (gjjtf + 9 P ^b(r) + -(1 + r 3 N )A(f) ) 7 o]^(r) = 0, (4) 

(-V r 2 + m>(f) = -[?¥jM] Pa tf = g„C N (a) Ps (r), (5) 

(-V 2 . + m 2 Juj(f) = g w p B (r), (6) 

(-V r 2 + m>(f) = |p 3 (r1, (7) 

(-V r 2 )A(f)=ep p (f), (8) 

where, p s (r), ps(f), Pz(f) and p p (r) are the scalar, baryon, third component of isovector, 
and proton densities at position f in the nucleus [20,21]. On the right hand side of Eq. (5), 
— [dMtf(a)/da] = g a CN{^)i where g a = g a (o = 0), is a new, and characteristic feature of 
QMC beyond QHD [28]. The effective mass for the nucleon, M^, is defined by: 

= ~n q gl J bag d 3 x ^(f)^(f) = -n q glS N (a) = - A [g a (a)a] , (9) 

where n q is the number of light quarks (u and d), and the MIT bag model quantities and 
the in-medium bag radius satisfying the mass stability condition are given by [19-21]: 

m* n (°) = E + U R ») 3B > ( 10 ) 

q=u,d U N 6 

s N (a) = [n* q /2 + m* q R* N (n; - 1)] / [n;(n; - 1) + m* q R* N /2] , (11) 

= y/tf + W"# 2 > m* q = m q - gla(r) , (12) 
dM* N /dR N \ RN=R% =0, (13) 

where g q a is the quark-a meson coupling constant. Here, the MIT bag model quantities are 
calculated in a local density approximation using the spin and spatial part of the wave func- 
tions, ip q (x) = N q e~ ttqt / R Nil) q {x) ) where N q is the normalization factor. The wave functions 
ip q (x) satisfy the Dirac equations for the quarks in the nucleon bag centered at position fin 
the nucleus (\x - fj < R* N [20,21]): 

<7 " d x - K - V?(fO) =F 7° (Vj(r) ± \v*tf))\ ( ) = , (14) 

where we approximate the constant, mean meson fields within the bag and neglect the 
Coulomb force. The constant, mean-field potentials within the bag centered at fare defined 
by V q (r) = gQa(f)i VJ(r) = g%uj(f) and V!*(r) = g q p b(r). The eigenenergies in units of l/R* N 
are given by: 



^Rn{v^±\v^). (15) 



In Eqs. (10) - (13), zn, B, x q , and m q are the parameters for the sum of the cm. and gluon 
fluctuation effects, bag pressure, lowest eigenvalues for the quark q, and the corresponding 
current quark masses, respectively. z N and B are fixed by fitting the nucleon mass in free 
space. We use the current quark masses m q=u> d = 5 MeV, and obtained = 3.295 and 
B = (170.0 MeV) 4 by choosing the bag radius for the nucleon in free space Rn = 0.8 fm. The 
parameters at the hadronic level, which are already fixed by the study of nuclear matter 
and finite nuclei [21], are as follows: m w = 783 MeV, m p = 770 MeV, m CT = 418 MeV, 
e 2 /47r = 1/137.036, g 2 a /An = 3.12, gl/An = 5.31 and gf/An = 6.93. The sign of m* in 
the nucleus in Eq. (12) reflects nothing but the strength of the attractive, negative scalar 
potential, and thus the naive interpretation of the mass for a physical particle, which is 
positive, should not be applied. 

At the hadronic level, the entire information on the quark dynamics is condensed into 
the effective coupling Cjv(er) of Eq. (5). Furthermore, when CV(c) = 1, which corresponds 
to a structureless nucleon, the equations of motion given by Eqs. (4)-(8) can be identified 
with those derived from QHD [28], except for the terms arising from the tensor coupling 
and the non-linear scalar field interaction introduced beyond naive QHD. 




As examples, we show in Fig. 1 the charge density distributions calculated for 40 Ca 
and 208 Pb, and also the energy spectra obtained for 40 Ca and 208 Pb in Figs. 2 and 3 [21], 
respectively. 

Next, we consider the limit of infinite symmetric nuclear matter [19-21]. In this limit all 
meson fields become constant, and we denote the mean- values of the u and a fields by ID and 
a. Then, equations for the cJ and self-consistency condition for the a are given by [19-21], 

"=tA-J d 3 k6(k F — k) — = ^p B , (16) 



(27r) 3 7 v ' ' ml37r 2 .- 2 



f d 3 k6(k F - k) Zm = = ^C N (a)p s , (17) 
( 27r ) J ^M* N 2 (a) + k 2 m l 



5 



where g„ = 3q^Sn(0) (see Eq. (11)), kp is the Fermi momentum, ps and p s are the baryon 
and scalar densities, respectively. Note that M^(cf) in Eq. (17) must be calculated self- 
consistently in the MIT bag model through Eqs. (9)-(14) for a given baryon density. This 
self-consistency equation for the a is the same as that in QHD, except that in the latter 
model one has Cn(p) = 1 [28]. 
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FIG. 2. Energy spectrum for 40 Ca [21] in the QMC model compared with experiment (Exp.), 
and that of QHD [28]. 
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FIG. 3. Same as Fig. 2 (for 209 Pb). 
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III. NUCLEAR MEDIUM MODIFICATION OF FORM FACTORS 



In this Section we outline the medium modification of the electromagnetic form fac- 
tors of the nucleon, as suggested in the recent polarization transfer measurements in the 
4 He(e, e'p) 3 H reaction [1,2]. The first data were analyzed in Ref. [1] using a variety of mod- 
els, nonrelativistic and relativistic, based on conventional nucleon-nucleon potentials and 
well-established bound state wave functions, including corrections from meson exchange 
currents, final state interaction and other effects [11-14]. The observed deviation, which 
was of order 10%, could only be explained by supplementing the conventional nuclear de- 
scription with the effects associated with the medium modification of the nucleon internal 
structure calculated by the QMC model [15-18]. 

In Fig. 4 we show the "super ratio", R/Rpwia, which was made for the final analysis of 
the polarization transfer measurements on 4 He [2]. Here, Rpwia stands for the prediction 
based on the relativistic plane- wave impulse approximation (PWIA) , and the measured ratio 
R is defined by: 



R 



(18) 
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FIG. 4. Super ratio R/ RpwiA, as a function of Q 2 , taken from Ref. [2]. See caption of Fig. 1 in 
Ref. [2] for detailed explanations. 



In Fig. 4, the modification of electromagnetic form factors of the bound nucleon cal- 
culated in the QMC model [15-18] (the solid line denoted by "Udias RDWIA + QMC") 
uses the improved cloudy bag model (ICBM) [29,30] for the free nucleon form factors. The 
ICBM [30] includes a Peierls-Thouless projection to account for center of mass and recoil 
corrections, and a Lorentz contraction of the internal quark wave functions. 

The electromagnetic current is given by the sum of the contributions from the quark core 
and the pion cloud, 

f( x ) = Q q e^ q {x)YM%) ~ ie[^{x)d^(x) - n(x)d^{x)} , (19) 
i 
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where Q q is the charge operator for a quark flavor q, and ir(x) destroys a negatively charged 
(or creates a positively charged) pion. Relevant diagrams included in the calculation of free 
space electromagnetic form factors are depicted in Fig. 5 [30]. 



N 2> N_ N B ^) C N N / \ N 

B 

(a) (b) (c) 

FIG. 5. Diagrams included in the calculation of free space electromagnetic form factors in the 
ICBM [30]. The intermediate baryon states B and C are restricted to the N and A. 



In the Breit frame the quark core contribution to the electromagnetic form factors of the 
bound nucleon is given by [15-18]: 

G* E (Q 2 )=V 2 Gt*(v 2 Q 2 ) , (20a) 
G* M (Q 2 )= V 2 Gt*(v 2 Q 2 ) , (20b) 

where Q 2 = -q 2 = q 2 , and the scaling factor r] = {M* N /E* N ), with E* N = ^ M* N 2 + Q 2 /A the 

energy and the mass of the nucleon in medium. G s ^ h ^ are the form factors calculated 
with the static spherical bag wave function, 

Gt*(Q 2 ) = ^Jd 3 r HQr) f q (r) K(r) , (21a) 
1 9 M* r 

Gt\Q 2 ) = p-^- j d 3 r UQt) (3* q j (x q r/R* N ) h{x q r/R* N ) K(r) . (21b) 

Here f q (r) = jl(x q r/R* N ) + (3* 2 j 2 (x q r/R* N ), where R* N is the nucleon bag radius in 
medium, x q the lowest eigenfrequency, and [3* 2 = (Q* — m*R* N )/(Q* + m*R* N ), with 

il* = \Jx 2 + {. m * q R*N) 2 and m* = m q — g q a u (see also Section II). The recoil function 
K{r) = jd 3 x f q {x)f q (—x — r) accounts for the correlation of the two spectator quarks, 
and D = jd 3 r f q (r) K{r) is the normalization factor. The scaling factor rj in the argument 
of m arises from the coordinate transformation of the struck quark, and the prefactor in 
Eqs. (20) comes from the reduction of the integral measure of the two spectator quarks in 
the Breit frame. 

The contribution from the pionic cloud is calculated along the lines of Ref. [15-18]. 
Although the pion mass would be slightly smaller in the medium than in free space, we use 
m* = m n , which is consistent with chiral expectations and phenomenological constraints. 
Furthermore, since the A isobar is treated on the same footing as the nucleon in the CBM, 
and because it contains three ground state light quarks, its mass should vary in a similar 
manner to that of the nucleon in the QMC model. As a first approximation we therefore take 
the in-medium and free space N-A mass splittings to be approximately equal, — ~ 
M A - M N . 
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The change in the ratio of the electric to magnetic form factors of the proton from free 
to bound, 

<„«vw)=(ig|)/(iis). (*> 

is illustrated in Fig. 6 for 4 He, 16 O (left panel) and for nuclear matter densities, p = po and 
P — \p$ (right panel) with po = 0.15 fm -3 . 




<2 2 (GeV 2 ) 

FIG. 6. The change in the ratio of the electric to magnetic form factors from free to bound 
protons, R^m/^em = (^e / G P m) / {G P E / G P M ) , for 4 He (lsi/ 2 state) and 16 (lsi/ 2 and ls 3 / 2 states) 
(left panel) [18] and nuclear matter (right panel) [25]. 16 0(B) stands for the results which allow 
changes of the bag constant according to the nuclear density (po = 0.15 fm -3 [17]). 

Because the average nuclear densities for all existing stable nuclei heavier than deuterium 
lie in the range ^po ~ P ~ Po, we consider these two specific nuclear densities to give the 
upper and lower bounds for the change of the electromagnetic form factors (and structure 
functions at large x) of the bound nucleon. We emphasize that in the present analysis the 
absolute value of the proton magnetic form factor at Q 2 = (the magnetic moment), which 
is enhanced in medium, plays an important role - as it did in the analysis of polarized (e, e'p) 
scattering experiments. 

Because of charge conservation, the value of G P E at Q 2 = remains unity for any p. On 
the other hand, the proton magnetic moment is enhanced in the nuclear medium, increasing 
with p, so that R P E * M < R P em a ^ Q 2 = 0- In fact, the electric to magnetic ratio is ~ 5% 
smaller in medium than in free space for p = |po, an d ~ 10% smaller for p = p . The effect 
increases with Q 2 out to ~ 2 GeV 2 , where the (bound/free) ratio deviates by ~ 20% from 
unity. 

The extension to the in-medium modification of the bound nucleon axial form factor 
G* A (Q 2 ) can be made in a straightforward manner [27]. Since the induced pseudoscalar 
form factor, Gp(Q 2 ), is dominated by the pion pole, and can be derived using the PCAC 
relation [29], we do not discuss it here. The relevant axial current operator is then simply 
given by 
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(23) 



where ip q (x) is the quark field operator for flavor q. 

Similarly to the case of electromagnetic form factors, in the preferred Breit frame the 
resulting bound nucleon axial form factor is given by [27]: 



g a (q 2 ) = v 2 gT*(v 2 Q% 

Gt\Q 2 ) = I U*r{\j%{x q r/R* N ) - i3fj!(x q r/R* N )]j (Qr) 



+2(3*y i (x q r/R* N )\j 1 (Qr)/Qr}} K{r)/D. 



(24) 



(25) 



In Fig. 7 we show the (normalized) free space axial form factor Ga{Q 2 ) calculated in the 
ICBM [27] together with the experimental data (left panel), and (the space component of) 
that calculated at nuclear densities p = (0.5, 0.7, 1.0, 1.5)p with p = 0.15 fm -3 . 




FIG. 7. Free space (normalized) axial form factor Ga{Q 2 ) calculated in the ICBM 
(left panel) [27] together with experimental data [31] summarized by a dipole form: 
Ga(Q 2 ) = 9a/ + Q 2 /m 2 A ) 2 with txia = (1-03 ± 0.04) GeV, and the ratio of in-medium to free 
axial form factors [27] (right panel), where qa = 1-14 is used in the ICBM calculation. 



At Q 2 = the space component G* A {Q 2 — 0) = g A is quenched [4] by about 10 % at 
normal nuclear matter density. The modification calculated here may correspond to the 
"model independent part" in meson exchange language, where the axial current attaches 
itself to one of the two nucleon legs, but not to the exchanged meson [4] . This is because the 
axial current operator in Eq. (23) is a one-body operator which operates on the quarks and 
pions belonging to a bound nucleon. The medium modification of the bound nucleon axial 
form factor G* A (Q 2 ) may be observed for instance in neutrino- nucleus scattering, similar to 
that observed in the "EMC-type" experiments, or in a similar experiment to the polarization 
transfer measurements performed on 4 He [1,2,9]. However, at present the experimental 
uncertainties seem to be too large to detect such medium effect directly. We should also 
note that the medium modification of the parity- violating F3 structure functions of a bound 
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nucleoli in deep-inelastic neutrino induced reactions can be extracted using the calculated 
in-medium axial form factor G* A (Q 2 ) and quark-hadron duality, as we discuss in the next 
Section for the case of electromagnetic form factors and the F 2 structure function. 

IV. QUARK-HADRON DUALITY AND NUCLEON STRUCTURE FUNCTIONS 

IN MEDIUM 

A new and interesting avenue for exploring medium modifications of hadronic observ- 
ables such as form factors and structure functions is provided by quark-hadron duality. The 
relationship between form factors and structure functions, or more generally between inclu- 
sive and exclusive processes, has been studied in a number of contexts over the years. Drell 
& Yan [32] and West [33] pointed out long time ago that, simply on the basis of scaling 
arguments, the asymptotic behavior of elastic electromagnetic form factors as Q 2 — > oo can 
be related to the x — > 1 behavior of deep-inelastic structure functions. In perturbative QCD 
language, this can be understood in terms of hard gluon exchange: deep-inelastic scattering 
at x ~ 1 probes a highly asymmetric configuration in the nucleon in which one of the quarks 
goes far off-shell after the exchange of at least two hard gluons in the initial state; elastic 
scattering, on the other hand, requires at least two gluons in the final state to redistribute 
the large Q 2 absorbed by the recoiling quark [34]. 

More generally, the relationship between resonance (transition) form factors and the 
deep-inelastic continuum has been studied in the framework of quark-hadron, or Bloom- 
Gilman, duality: the equivalence of the averaged structure function in the resonance region 
and the scaling function which describes high W data. The recent high precision Jefferson 
Lab data [23] on the F 2 structure function suggests that the resonance-scaling duality also 
exists locally, for each of the low-lying resonances, including surprisingly the elastic [24], to 
rather low values of Q 2 . 

A number of recent studies have attempted to identify the dynamical origin of Bloom- 
Gilman duality using simple models of QCD [35-37]. It was shown, for instance, that in a 
harmonic oscillator basis one can explicitly construct a smooth, scaling structure function 
from a set of infinitely narrow resonances [35,36]. Whatever the ultimate microscopic origin 
of Bloom-Gilman duality, for our purposes it is sufficient to note the empirical fact that 
local duality is realized in lepton-proton scattering down to Q 2 ~ 0.5 GeV 2 at the 10-20% 
level for the lowest moments of the F 2 structure function. In other words, here we are not 
concerned about why duality works, but rather to assess the phenomenological consequences 
of the fact that it does work. 

Motivated by the experimental verification of local duality, one can use measured struc- 
ture functions in the resonance region to directly extract elastic form factors [38] . Conversely, 
empirical electromagnetic form factors at large Q 2 can be used to predict the x — > 1 behav- 
ior of deep-inelastic structure functions [22,39,40]. The assumption of local duality for the 
elastic case implies that the area under the elastic peak at a given Q 2 is equivalent to the 
area under the scaling function, at much larger Q 2 , when integrated from the pion threshold 
to the elastic point [22]. Using the local duality hypothesis, de Rujula et al. [38], and more 
recently Ent et al. [24], extracted the proton's magnetic form factor from resonance data 
on the F 2 structure function at large x, finding agreement to better than 30% over a large 
range of Q 2 (0.5 < Q 2 < 5 GeV 2 ). In the region Q 2 ~ 1-2 GeV 2 the agreement was at the 
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~ 10% level. In Fig. 8 we show the extracted proton magnetic form factor G P M using the 
quark-hadron local duality relation [41] and the F p (£) parametrization of Ref. [42], in order 
to estimate how reliable the quark-hadron duality assumption may be. 




Q 2 (GeV 2 ) 

FIG. 8. Duality extraction of the free space proton magnetic form factor [41] using the 
quark-hadron duality relation and the parametrization in Ref. [42]. 



Applying the argument in reverse, one can formally differentiate the local elastic duality 
relation [22] with respect to Q 2 to express the scaling functions, evaluated at threshold, 
x = x t h = Q 2 / (W 2 h - Mjj + Q 2 ), with W t h = M N + m n , in terms of Q 2 derivatives of elastic 
form factors. In Refs. [22,39] the x — > 1 behavior of the neutron to proton structure function 
ratio was extracted from data on the elastic electromagnetic form factors. Extending this 
to the case of bound nucleons, one finds that as Q 2 — > oo the ratio of bound to free proton 
structure functions is: 

FT dGV 2 /dQ 2 

F$ dG^ldQ 2 ' 1 ' 

At finite Q 2 there are corrections to Eq. (26) arising from G P E and its derivatives, as discussed 
in Ref. [39]. (In this analysis we use the full, Q 2 dependent expressions [39,40].) Note that 
in the nuclear medium, the value of x at which the pion threshold arises is shifted: 

* _ ( m n {2M N + m^) + Q 2 \ 
Xth " Xth " {m«(2(M* N + V N )+m n ) + Q 2 ) Xth > {Z<) 

where Vn = 3^ uj is the vector potential felt by the nucleon and (consistent with chiral 
expectations and phenomenological constraints) we have set m* = m n . (See also Eq. (16).) 
However, the difference between x t h and x* th has a negligible effect on the results for most 
values of x considered. 

Using the duality relations between electromagnetic form factors and structure functions, 
in Fig. 9 we plot the ratio F^/F^ as a function of x, with x evaluated at threshold, x = Xth 
(solid lines). 
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FIG. 9. The ratios of the bound to free proton structure functions F%* /F%, calculated in the 
QMC and PLC models. 



We emphasize that since we are interested in ratios of form factors and structure functions 
only, what is more relevant for our analysis is not the degree to which local duality holds for 
the absolute structure functions, but rather the relative change in the duality approximation 
between free and bound protons. Note that at threshold the range of Q 2 spanned between 
x = 0.5 and x = 0.8 is Q 2 « 0.3-1.1 GeV 2 . Over the range 0.5 < x < 0.75 the effect is almost 
negligible, with the deviation of the ratio from unity being p$ 1% for p = |po and ^ 2% 
for p = po- For x ^ 0.8 the effect increases to ~ 5%, although, since larger x corresponds 
to larger Q 2 , the analysis in terms of the QMC model is less reliable here. However, in the 
region where the analysis can be considered reliable, the results based on the bound nucleon 
form factors inferred from the polarization transfer data [1,2] and local duality imply that 
the nucleon structure function undergoes very little modification in medium. 

It is instructive to contrast this result with models of the EMC effect in which there 
is a large medium modification of nucleon structure. For example, let us consider the 
model of Ref. [26], where it is assumed that for large x the dominant contribution to the 
structure function is given by the point-like configurations (PLC) of partons which interact 
weakly with the other nucleons. The suppression of this component in a bound nucleon is 
assumed to be the main source of the EMC effect. This model represents one of the extreme 
possibilities that the EMC effect is solely the result of deformation of the wave function 
of bound nucleons, without attributing any contribution to nuclear pions or other effects 
associated with nuclear binding [44] . 

The deformation of the bound nucleon structure function in the PLC suppression model 
is governed by the function [26]: 

5(k) = l-2(k 2 /2M + e A )/AE A , (28) 

where k is the bound nucleon momentum, e A is the nuclear binding energy, and AE A ~ 0.3- 
0.6 GeV is a nucleon excitation energy in the nucleus. For x ^ 0.6 the ratio of bound to 
free nucleon structure functions is then given by [26]: 
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(29) 



The x dependence of the suppression effect is based on the assumption that the point-like 
configuration contribution in the nucleon wave function is negligible at x ^ 0.3 (FJ?*/FJ? = 
1), and for 0.3 ^ x ^ 0.6 one linearly interpolates between these values [26]. The results for 
4 He and 16 are shown in Fig. 9 (dashed lines) for the average values of nucleon momentum, 
(k 2 ), in each nucleus. The effect is a suppression of order 20 % in the ratio F?*/F? for 
x ~ 0.6-0.7. In contrast, the ratios extracted on the basis of duality, using the QMC model 
constrained by the 4 He polarization transfer data [1,2], show almost no suppression (< 1- 
2%) in this region. Thus, for 4 He, the effect in the PLC suppression model is an order of 
magnitude too large at x ~ 0.6, and has the opposite sign for x > 0.65. 

Although the results extracted from the polarization transfer measurements [1,2] rely on 
the assumption of local duality, we stress that the corrections to duality have been found to 
be typically less than 20% for 0.5 < Q 2 < 2 GeV 2 [23,43]. The results therefore appear to 
rule out large bound structure function modifications, such as those assumed in the point- 
like configuration suppression model [26] , and instead point to a small medium modification 
of the intrinsic nucleon structure, which is complemented by standard many-body nuclear 
effects. 

Nevertheless, given the large differences between the theoretical predictions, data on 
nuclear structure functions at large x would be extremely valuable in discriminating between 
these scenarios. Possible insights into the medium modifications of bound nucleon structure 
functions may be garnered from semi-inclusive deep-inelastic scattering experiments from 
deuteron, by studying the spectrum of tagged spectator nucleons [45]. As a consistency check 
on the analysis, one can also examine the change in the form factor of a bound nucleon that 
would be implied by the corresponding change in the structure function in medium. Namely, 
from the local duality relation [38,40]: 

[ G " (Q2) ^^(^£™- (30) 

one can extract the magnetic form factor by integrating the -Ff (£) structure function over 
£ between threshold, £ = £ t h, and £ = 1. Here £o = C( x — 1)> A*p is the proton magnetic 
moment, and £ = 2x/(l + \Jl + x 2 /r) with r = Q 2 /4M^ and x the Bjorken variable. 
In Fig. 10 we show the PLC model predictions for the ratio of the magnetic form factor 
of a proton bound in 4 He to that in vacuum, derived from Eqs. (29) and (30), using the 
parameterization for F$(£) from Ref. [24], and an estimate for the in-medium value of n* from 
Ref. [15-18]. Taking the average nucleon momentum in the 4 He nucleus, k — (k) « 135 MeV, 
the result is a suppression of about 20% in the ratio G P ^/G P M at Q 2 ~ 1-2 GeV 2 (solid 
curve). Since the structure function suppression in the PLC model depends on the nucleon 
momentum (Eq. (28)), we also show the resulting form factor ratio for a momentum typical 
in the (e,e'p) experiment, k = 50 MeV (long dashed). As expected, the effect is reduced, 
however, it is still of the order 15% since the suppression also depends on the binding energy, 
as well as on the nucleon mass, which changes with density rather than with momentum. 
In contrast, the QMC calculation, which is consistent with the MAMI 4 He quasi-elastic 
data, and Jlab polarization transfer measurements on 4 He [1,2], produces a ratio which is 
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typically 5-10% larger than unity (dashed). Without a very large compensating change 
in the in-medium electric form factor of the proton (which seems to be excluded by y- 
scaling constraints), the behavior of the magnetic form factor implied by the "PLC model + 
duality" would produce a large enhancement of the polarization transfer ratio, rather than 
the observed small suppression [1,2] (see Eq. (1) and Fig. 4). 



1.2 




0.7 1 ■ ' ■ 1 ■ ' ■ 1 

0.5 1 1.5 2 

G 2 (GeV 2 ) 

FIG. 10. The ratios of the bound proton to free magnetic form factors G v ^jG v M ^ calculated in 
the QMC and PLC models. 



V. SUMMARY 

In this article we have discussed the medium modification of the internal structure of 
bound nucleons due to the change of quark response to the nuclear environment. Recent 
experimental results [1,2] from polarized proton knockout reactions off 4 He nuclei have pro- 
vided the first direct evidence for a possible small but nonzero modification of the proton 
electromagnetic form factors in the nuclear medium. The analyses in Refs. [1,2] found that 
when compared with conventional nuclear calculations, the medium modifications observed 
in the 4 He data could only be accounted for within models in which the conventional de- 
scriptions are supplemented by a small modification of the nucleon form factors in medium 
calculated in the QMC model [15-21] (see also Ref. [46]). 

We have also analyzed the medium modification of the bound nucleon axial form factor 
Ga(Q 2 )- The modification of the axial form factor may be observed for instance in future 
high precision neutrino-nucleus scattering experiments [47]. 

Finally, we have examined the consequences of quark-hadron duality applied to nucleons 
in the nuclear medium. Utilizing the experimental results [1,2] we have used local duality to 
relate model-independently the medium modified electromagnetic form factors to the change 
in the intrinsic structure function of a bound proton. While the results rely on the validity of 
quark-hadron duality, the empirical evidence suggests that for low moments of the proton's 
F 2 structure function the duality violations due to higher twist corrections are < 20% for 
Q 2 > 0.5 GeV 2 [23], and decrease with increasing Q 2 . In the context of the QMC model, the 
change in nucleon form factors allowed by the data imply a modification of the in-medium 
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structure function of ^ 1-2% at 0.5 ^ x < 0.75 for all nuclear densities between nuclear 
matter densities, p = po, and p = |p - 

The results place rather strong constraints on models of the nuclear EMC effect, espe- 
cially on models which assume that the EMC effect arises from a large deformation of the 
nucleon structure in medium. While suggesting the need for explicit quark degrees of free- 
dom in the nucleus, our findings appear to disfavor models with large medium modifications 
of structure functions as viable explanations for the nuclear EMC effect, although it would 
be desirable to have more data on a variety of nuclei and in different kinematical regions. 
A such proposed experiment [48] at Jefferson Lab on 16 at Q 2 = 0.8 GeV 2 , which would 
make use of other, high precision cross section data at this momentum transfer, would have 
about 15 times the statistics of the original commissioning experiment [9]. This would en- 
able a more thorough comparison of the medium dependence of form factors and structure 
functions for different nuclei. 
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